Variation of dispersion measure: evidence of geodetic precession of binary pulsars 
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Variations of dispersion measure (DM) have been observed in some binary pulsars, which can 
not be well explained by the propagation effects, such as turbulence of the interstellar media (ISM) 
between the Earth and the pulsar. This paper provides an alternative interpretation of the phe- 
nomena, the geodetic precession of the orbit plane of a binary pulsar system. The dynamic model 
can naturally avoid the difficulties of propagation explanations. Moreover the additional time delay 
represented by the DM variation of two binary pulsars can be fitted numerically, through which 
some interesting parameters of the binary pulsar system, i.e., the moment of inertia of pulsars can 
be obtained, Ii = (2.0 ± 0.6) x 10 45 g cm 2 . The elimination of the additional time delay by the 
dynamic effect means that ISM between the these pulsars and the Earth might also be stable, like 
some other binary pulsars. 



PACS numbers: 



I. INTRODUCTION 



DM within ISM can delay a radio pulse in reaching 
Earth by a number of seconds equal to DM/(2A\ x 
10~ 4 / 2 ), where / is the observing frequency in MHz and 
DM is the column density of free electrons integrated 
along the line of sight in unite of pc cm -3 [[J, 



DM 



ripdz . 



(1) 



where I is the distance to the pulsar. For many pulsars, 
the DM can be characterized as a constant that holds 
steady over years of observation. 

However, millisecond pulsar in the globular cluster 47 
Tucanae, i.e., PSR J0023-7203J (47 Tuc JVshows vari- 
ations of DM as a function of orbital phase 2] . The vari- 
ations of DM are independent of frequency, which indi- 
cates that the additional time delay is not likely a propa- 
gation effect 0, since propagation effect predicts that the 
waves at the low-frequency and high-frequency should 
show very different time of arrivals (TO As). 

The galactic binary pulsar PSR J0621+1002 experi- 
ences dramatic variability in its DmQ , with gradients 
as steep as 0.013 pc cm" 3 yr -1 . If the DM variation is 
interpreted as spatial fluctuation in the interstellar elec- 
tron density, then it would obviously deviate from the 
simple power law predicted by the standard theories of 
ISM[!||. 

Therefore, as discussed by the authors 0, Q, attribut- 
ing the additional time delay (or residuals of (TOAs) to 
propagation effect, DM variation, is not very satisfactory 
in the comparison with the observations. 

The geodetic precession induced orbital effect of a bi- 
nary pulsar system can cause an additional time delay, 
which can well explain the long-term variabilities, such 
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as derivatives of the semi- major axis, x, x, and the or- 
bital period, P b P b measured in PSR J2051-0827 and 
PSR B1957+2O0. 

This paper applies the geodetic precession induced 
time delay in shorter time scales (relative to secular vari- 
abilities) to interpret the residuals in timing measure- 
ment of 47 Tuc J, which has been attributed to the vari- 
ation of DM. The new explanation can fit the residuals 
and also avoid the frequency difficulty in the propagation 
effect. 

Moreover, the geodetic precession induced variations 
at different time scales impose strong constraints on the 
intrinsic parameters of a binary pulsar system. Fitting 
them together, we can obtain for the first time numerical 
result of the spin angular momenta of the two stars as well 
as the moment of inertia of the pulsar which is consistent 
with theoretical predictions. 

The significant DM variation of PSR J0621+1002 and 
PSR J0024-7204H (47 Tuc H) can also be well explained 
by the dynamic effect. The elimination of the additional 
time delay by the dynamic effect means that the DM (or 
ISM) of these three binary pulsars might be very stable. 

The situation is similar to PSR B1855+09, which has 
found no unexplained perturbation in the long record 
of TOAs, and lead to a reduction in the upper limit 
of energy density in the gravitational wave background 
radiation 0. 

In section II the geodetic precession induced orbital 
precession velocity in general case is introduced. In sec- 
tion III the additional time delay due to geodetic preces- 
sion of a binary pulsar system is derived. And in section 
IV, V and VI the dynamic effect is applied to 47 Tuc 
J, PSR J0621+1002 and 47 Tuc H respectively. Section 
VII summarizes the relation of the geodetic precession 
induced effects with the time delay in two cases, and also 
the evidences the geodetic precession in binary pulsars. 
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II. ORBITAL PRECESSION 

The motion of a binary system can be regarded as the 
precession of three vectors, the spin angular momenta 
of the pulsar and its companion star, Si and S2, and 
the orbital angular momentum L. The change of the 
orbital period due to the gravitational radiation is 2.5 
post-Newtonian order (2.5 PPN), whereas the geodetic 
precession corresponds to 1.5 PPN. So the influence of 
gravitational radiation on the motion of a binary sys- 
tem can be ignored when discussing dynamics of a binary 
pulsar system. Therefore, the total angular momentum, 
J = L + Si + S2, can be treated as invariable both in 
magnitude and direction (J = 0). With Slo denoting the 
precession rate of L_around J, the spin-orbit coupling can 
be expressed as[~ 



f2o x L = — Jli x Si — tt 2 x S 2 



(2) 



where fii and fi 2 represent the precession of the pulsar 
and its companion star, respectively. Ignoring terms over 
2 PPN, fii and fi 2 can be written as[7| 



2r° mi 



^(4+^), (3) 



2r 



m 2 



where mi and m 2 are masses of the pulsar and the com- 
panion star respectively, and r is the separation of mi 
and m 2 . Notice L ~ r 1/2 , fii and fi 2 are 1.5 PPN. 

Barker and O'Connell's two-body equation included 
two spins, but the orbital precession velocity was not 
expressed as relative to the total angular momentum, J, 
therefore, it cannot be compared to observation directly 
(J is static relative to the line of sight after counting out 
the proper motion of the binary system). 

Apostolatos et al and Kidder [8J, |2j's orbital precession 
velocity was relative to J, however their velocity of orbit 
plane was derived in the case of one spin, i.e., Si = 0. 
Which is suitable only for special binary systems, like 
pulsar-black hole binary. 

Therefore, it seems contradictory that in Barker and 
O'Connell's equation L doesn't precess around J, but 
in practical use, L is expressed as precessing around 
J. Actually these two expressions can be consistent in 
the scenario which has been mentioned by Smarr and 
BlandfordllO|. In which L, Si and S 2 all precess around 
J rapidly (1.5PPN), whereas the velocities of L relative to 
Si and S 2 are very small (2PPN). And in the confronta- 
tion with observation, only the rapid precession velocity 
relative to J, 1.5PPN, should be used. 

GongQ derived the orbital precession velocity in a 
general cases, which is relative to J, and includes both 
Si and S 2 . The derivation is based on two simple as- 
sumptions: conservation of the total angular momen- 
tum, EqQ, and the geometry constraints of the triangle 
formed by J = S+L (S = Si + S 2 ). The two assumptions 
lead to precession rate of L around J [f| , 



fin 



n 2 sinAis + (fii 



S 11 

^2)^ sinALSi 



(4) 



where 5}' = Si cos rjssi , denoting the component of Si in 
the plane determined by S and J. Note that LsinA^j w 
S is used in Eq|0J, since S/L <C 1. The right-hand side 
of EqQJ can as well be written by replacing subscribes 1 
with 2 and 2 with 1. 

The geodetic precession of the orbit can cause an ad- 
ditional apsidal motion. In the case of S/L <C 1, the 
advance of the precession of the periastron, Co can be 
given by [T(l | 



• obs 

LO = LO 



GR 



fi COS X L J 



LO GR 



fin 



(5) 



As given by Eq|@}, fi can be as large as Co GR (1.5 PPN), 
the GR prediction of the advance of periastron. 

The effect of spin-orbit coupling on secular evolution 
of the orbital inclination, i, can be given by 



cos A lj cos / — sin Al j sin I cos r/n 



(6) 



where I is the angle between the total angular momen- 
tum, J, and the line of sight, and t/q = £l t + rji (m is the 
initial phase) is the phase of precession of L. Thus i is 
also a function of time. By EqlJSJ the first derivative of 
the projected semi-major axis is[5j 



= — .xfin sin Xlj sin rj cot i 



(7) 



However, since Si and S2 precess with different veloci- 
ties, fii and fi 2 respectively (mi ^ m 2 ), then S varies 
in both magnitude and direction (Si, S2 and S form a 
triangle), then from the triangle of S, L and J, in react to 
the variation of S, L must vary in direction (|L| =const), 
which means the variation of Xlj (J is invariable). 

The change of Xlj means that the orbital plane tilts 
back and forth. In turn, both Xls and Xjs vary with 
time. Therefore, from Eq(0}, the derivative of the rate of 
orbital precession can be given byQ, 

fi = fi 2 fii 2 A 3 A4 — fii 2 Ai(fi iA 2 + fii 2 A 3 ) , (8) 



fii - fi 2 , fim = fii - fi , Ai = -J- sin X LSl , 
X 



= S V , S V2 Bin.M.2 and X = cos- A 



where fii 2 
A 2 = tan 77^1 

with a = sin Xjs + ^rff ' Svi = Sl sin Xjsi and Sv2 = 
S2 sin A jsi represent components of Si and S2 that are 
vertical to J. 

Note that fii and fi 2 are unchanged when ignoring the 
orbital decay (2.5PPN), and Als q are unchanged also, 
since they decay much slower than the orbital decay 
fin can be easily obtained through EqJHJ). 

fin, the derivative of fio, can be absorbed by iV The 
variation in the precession velocity of the orbit results in 
a variation of orbital frequency (vb — 2iv/Pb), v' h — v\, = 
fioAi . Then we have i>b = fin, therefore 



Pb = - 



2tt 



(9) 



From Eq© and Eq@, we can see that the contribution 
of fin to Pb can be as large as 1 PPN, which is much larger 
than the contribution of GR to P b (2.5PPN). While P b 
can also be much smaller than 1 PPN in special combi- 
nation of parameters in Eq© . 
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III. GEODETIC PRECESSION INDUCED TIME 
DELAY 



As discussed in section II, the geodetic precession in- 
duced orbital effect results an additional apsidal motion 
which can be absorbed by the post-Kepler parameter, 
, an additional precession of orbital plane which can 



,obs 



be absorbed by, i, and an additional variation of orbital 
period which can be absorbed by PbQ- 

These additional effects of a binary system can not only 
cause long-term (secular) time delay, but also short-term 
time delay. 

The essential transformation relating solar system 
barycentric time tf, to pulsar proper time T is summa- 
rized by the expression 



t b - t Q = T + A R + A E + A s + A A , 



(10) 



where Ar is the " Roemer time delay" , is the propagation 
time across the binary orbit; and Ag are the orbital 
Einstein and Shapiro delays; and Aa is a time delay re- 
lated with aberration caused by rotation of the pulsar. 
The dominant time delay, Ar is given 01 



where 



A R = xF(co + u) , 



F(io + u) = sinu;[cosM — e(l + S r )] + 



(11) 



[l-e 2 (l-fc) a ] 1/2 cosw 



sin u . 



(12) 



In calculation the small quantities, 5 r and Sg due to aber- 
ration are ignored, u and A e {u) are the eccentric anomaly 
and the true anomaly respectively. The relations of it, 
A e (u) and the longitude of periastron, to are given |12j 



o Pb,T~T a 2 
u — esmu = 2ir[( ) (- 



Pb 



Pb 



(13) 



A e (u) = 2arctan[(i±£) 1 / 2 tan~], (14) 
1 — e 2 



to = luq + kA e (u) . 



(15) 



where k = CoPb/(2ir). Ea lfTU|) to Ea ((TB|l are the standard 
treatment in pulsar timing measurement. When consider 
the contribution of orbital precession to the time of ar- 
rival, the Roemer delay should has a different value, A^, 
relative to the standard, Ar, which doesn't include the 
dynamic effect. The deviation leads to an additional time 
delay. 

SA R = A' R — Ar = x'F{uj' + u') - xF(uj + u) . (16) 

x' = x + xt at right-hand side of Eg ||16|) corresponds to 
an additional time delay by the precession of the orbital 



plane, in which x is given by EqJTJ. u' represents an 
additional time delay by the nutation, in which of 
Eci H13Jl ) should be given by Eq©). And oj' represents 
an additional time delay by the apsidal motion, in which 
k of Ea ljT5|l should be written as kl — (Co + fl ) Pb / (2ir) . 

Thus SAr of Ea (1 16(1 at two different moment, t and 
t + t can be given as 



SA^ 



SA\ 



[x'F(oj' + u') - xF(oj + u)] 



t + T 



-\x'F(u)' + v!) -xF(u+u)] t 



(17) 



In other words if the true Roemer delay is given by A' R , 
which include the dynamic effect, but it is treated as the 
standard one Ar without the dynamic effect, then ad- 
ditional time delay results. Following section shows that 
the additional time delay can be well fitted by the differ- 
ence between A^, and A^ given by Ea (|16fl and Ea i|17|) . 



IV. 47 TUC J 



A. estimations 

This section shows that in the geodetic precession 
model, once an orbital precession velocity, Oo, which is 
suitable to explain one secular variability, i.e., x, then it 
is also suitable to explain all the other secular variabili- 
ties, like Co and Pb and DM of a binary pulsar. Which 
indicate that the physics underlying these phenomena is 
most likely the geodetic precession. 

With mi = 1.44M , m 2 = O.O2M , P b = 0.12d 
and x obs — 0.040, we have the semi-major axis, a — 
[GM /vl] 1 / 3 = 8.2 x 10 10 cm, and orbital angular momen- 
tum, L = v b [ia 2 (\ - e 2 ) 1/2 = 1.6 x 10 50 g cm 2 s _1 , with /x 
the reduced mass. Then we have f2i = 4.4 x 10 _11 s _1 , 
fi 2 = 2.5 x 10 _9 s _1 by Eq©. Assume Q « Qi, and by 
Ea 11611 the geodetic precession induced time delay in the 
time interval r = lyr is approximately, 

\5A R I w x\S sin(w)| ks xQ t\ cos(lu)\ w 5.5 x 10 _7 (s) , 

while in the case £Iq ~ O2, the additional time delay 
in one year is \Ar\ w 3.2 x 10~ 5 . The observed one 
is DM < 2 x 10 _4 cm~ 3 pc/yr|2|, which corresponds to 
2 x 10~ 6 s per yr. 

Therefore, if £7i < Oq < ^2, i.e., ~ ^2/10, then the 
corresponding time delay, |A^| « 3.2 x 10~ 6 s, is close to 
the observational limit. 

The geodetic precession induced secular variabilities, 
x and Pb can also be estimated. With £lo ~ ^i and 
± ob s = (__ 2 y ± 0.7) x 10~ 14 2], EqQ) becomes 



1.8 x 10 _12 sinA L j|sin77 cotz| w 2.7 x 10" 



(19) 



Thus sinALj| sinryocotil w 1.5 x 10 2 . If | sin 170 cot i| w 
0.1, then Xlj sa 1.5 x 10 — x , which means S/L w 1.5 x 
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10~ 4 . With L obtained above, we have Si ~ S ~ 2.4 x 
10 49 g cm 2 s -1 . Having the measured pulsar period, P = 
2.1ms, the moment of inertia of the pulsar is Ii sa 8.O/45 
(745 = 1 x 10 45 g cm 2 ). While if assuming £1 r* Q 2 , then 
the moment of inertia is I\ ~ 0.1 4/45. And similarly in 
the case fio ~ ^2/10, I\ « I.4/45. Which means the 
velocity that suitable for DM corresponds to a moment 
of inertia that is very close to the theoretical prediction. 

As shown in Eq@, P& is determined by f2o, which can 
vary in a range of several order of magnitude by different 
combination of variables in EqJSJ). If we assume Qq w 
n^l.9x lCT 21 s- 2 , then by Eq©, we have P b w 0.33 x 
10~ 13 . Which is about one order of magnitude smaller 
than the measured one. P§ bs = (-0.5 ± 0.13) x 10~ 12 . 
Whereas, if we assume f2 ~ fi 2 ., then Pb « 9.8 x 10 -11 . 
Which is two order of magnitude larger than P b ofcs . 

And in the case n ~ fi 2 /10, Pb ~ 9.8 x 10" 13 , which 
is close to the measured P b obs . 

Therefore, the orbital velocity of order of magnitude, 
f2o ~ ^2/10, can consistent with measured variabilities. 
Which indicates that the geodetic precession might the 
true mechanism that responsible for the observational re- 
sults. 



DM and A obs can be given byQ 



SA obs = 



DM(t 



DM(t) 



2.41 x 10" 4 / 2 2.41 x 10" 4 / 2 



(20) 



Thus the geodetic precession induced time delay at two 
moment t and t + r given by Ea i|17fl . can be used to ex- 
plain the observational one represented by Ea l|20|l , <5 A obs . 
The new interpretation can not only solve the difficulties 
in the previous explanations of 47 Tuc J, but also fit the 
additional time delay numerically. 

The estimation above indicates that SAp, x and Pb can 
be well explained by the geodetic precession of the binary 
system. Now we can fit these variations numerically. 

The vectors Si, S2 and S are studied in the coordi- 
nate system of the total angular momentum, in which 
the z-axis directs to J, and the x- and y-axes are in the 
invariance plane. S can be represented by Sp and Sy, 
the components parallel and vertical to the z-axis, re- 
spectively: 



S = (S v 



Sp) 1 ' 2 



(21) 



Sp and Sy can be expressed (recall Svi, Sv2 and Sy 
form a triangle) as 



Sp = Si cos Xjsi + S 2 cos Xjs2 



B. fitting 

DM variation with the orbital phase is clearly detected 
in 47 Tuc J, as shown by the scattered points with error 
bars in Fig 2. The variation has been interpreted as a 
cometary-like phenomenon, which caused by material at 
a considerable distance from the companion, with much 
higher average electron density than that of the rest area. 
In other words, the plasma cloud is responsible for the 
additional time delay, which can be attributed to the 
variations of DM . 

Whereas the variation of DM and the residuals of 
TOAs as function of orbital phase are very close at dif- 
ferent frequencies, from 660MHz to 1486MHz0. Which 
means that the additional time delay is not likely caused 
by the propagation effect of ISM. Since the time delay 
due to DM variations should be very different at low and 
high-frequency. 

Actually, what measured in 47 Tuc J of Fig 2 is resid- 
ual of TOAs, or additional time delay, A obs , which can 
be attributed to the variation of DM. And the discus- 
sion above indicates that such explanation is difficult to 
explain the frequency problem. 

Therefore, we can transform the DM variation of Fig 2 
back into the additional time delay, A ohs , and try to in- 
terpret it by the dynamic effect. 

The time delay, A obs , measured in the time interval 
between two moment, t and t + r (or orbital phases <p(t) 
and (j>{t + t)), corresponds to the DM variation, DM{t) 
and DM(t + t) respectively. The relationship between 



Sy = {Syi + S V2 - 2S V1 Sy 2 COS7 1S1S2 ) 1 / 2 , (22) 

where r]siS2 is the misalignment angle between Syi and 
Sy 2l which can be written as 

VS1S2 = (fti - n 2 )t + <Pt . (23) 

Therefore, by the variation of S as function of time (in 
the case of one spin, S =const), we can obtain Qq as 
function of time through EqQ . Thus the measured A ofcs 
of Ea lt2U|) (or DM) can be fitted step by step through 
Eqjljj. 

SA obs = (5A^ +T - SAp . (24) 

Through Eq© and Eq(J7| , the measured secular variabil- 
ities, P b obs and i obs 0, can also be fitted along with A ofcs 
by Monte-Carlo method. Obviously the long-term and 
short-term together imposes very stringent constraints on 
the numerical solutions. Notice that the measured mass 
function, fi(M Q ) = 4.864 x 10~ 6 @ is also considered in 
the fitting. 

As shown in Eq© and EqJ7J), fio is included in SAp 
and x, meanwhile f2o is included in Both Qq and 
f2o contain S(t), Si and angles, as shown in Eq(@J and 
EqlJSJ. Therefore, fitting the short-term and long-term 
parameters lead to the determination of Si and 62, and 
in turn the moment of inertia of the pulsar, Ii, since the 
pulsar period is known. 

In numerical fitting, Si and S2 are fitted in the range 
[2.5xl0 46 ,1.4xl0 49 ] and [2.5xl0 46 , 7.6xl0 49 ] (gcnrV 1 ) 
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respectively, which are enough to cover the estimated val- 
ues. The best solution is shown in Table I. By Eg (|2UI> and 
Eq|@J, we have DM oc 5A R oc fi oc Si, thus the 30% 
errors in the DM variation in the upper plot of Fig 2, 
can cause about 30% error in the fitted results, such as 
Si and I\. 

After fitting the measured time delay by the predicted 
one, as Eal|24fl. we can transform the predicted one, at 
the right hand side of Ea l|24l) into the theoretical DM 
variation, as displayed by the solid curves of Fig 2. So 
that it can be compared with the measured ones more 
clearly. 

The measured DM variation are slightly different at 
different time, i.e., 1998 June and 1999 October, as shown 
in Fig 2 (points with error bars), which can be explained 
by Ea l|16l) from which 5A R are slightly different at dif- 
ferent time. 

V. PSR J0621+1002 

For PSR J0621+1002, dramatic variability of its DM 
has been measured, with gradient as steep as 0.013 pc 
cm -3 yr -1 , shown in Fig 3. By the standard picture, the 
turbulence spreads energy from longer to shorter length 
scales arises a power law of the structure function, which 
is given by D DM (r) = {[DM(t + r) - DM(t)] 2 ), where r 
is the time lag between DM measurement. However the 
structure function obtained from observation obviously 
deviates from the simple power lawQ. Moreover, there 
is also no obvious differences in DM variation correspond- 
ing to 430 and 1410 MHz@. Which indicates that DM 
variation of this binary pulsar is also independent of fre- 
quency. In the geodetic precession induced model these 
difficulties can be explained naturally. 

Similarly, with mi = 1.70Mp, m 2 = 0.97M©, e = 
0.00246, P b = 8.3d and x = 120, we have fJ x = 1.4 x 
10- 12 s -1 , fl 2 = 2.2 x 10" 12 s _1 . Assuming fl w ft 2 , the 
geodetic precession induced time delay given by Ea l(TH|) 
can be written as (r = lyr) 

\SA R \ w xQ t\ cos(w)| w 8.3 x 10~ 4 (s) , (25) 

The measured maximum DM variation, DM — 0.013pc 
cm~ 3 yr _1 Q, corresponds to a delay of A ohs = 2.9xl0~ 4 s 
per year. Therefore the predicted one and the measured 
one can be consistent. 

VI. 47 TUC H 

The DM variation of 47 Tuc H, DM = 
-0.024(3)cm- 3 pc yr" 1 (equivalent A R = 5.4 x 10~ 4 s 
per year), is one of the largest DM variations for any 
pulsar, and no binary parameters available can explain 
the trend0- Whereas the geodetic precession model can 
explain this large DM variation easily. 

By the same treatment as the two binaries above, with 
mi = 1.50M o , m 2 = 0.186M Q and x obs = 2.150, the 



additional time delay per year is approximately, A R w 
1.8 x 10" 4 s, in the case fi ~ ^i ~ 2.7 x 10 -12 s -1 ; and 
A R « 1.2 x 10~ 3 s, in the case fl ~ ^2 ~ 1.7 x 10 -11 s -1 . 
Therefore, the large DM measured in 47 Tuc H, which 
is A R = 5.4 x 10~ 4 s per year, can be well explained 
by an orbital precession velocity that is in the range, 

Hi <n < ^2- 

VII. DISCUSSION 

As shown in Fig 4, the the orbital precession velocity, 

00, can contribute to the orbital precession, absorbed by 
X] apsidal motion, absorbed by u>; and nutation, absorbed 
by Pb respectively. 

For very special NS-NS or NS-WD binary pulsars, 
mi = m 2 or one spin is ignorable (i.e, Si = 0), S is 
a constant vector, then f^o is a constant, which means 
= 0, and in turn Pi, = 0, therefore, there will be only 
static orbital precession and apsidal motion but no nuta- 
tion. In such special cases the additional time delay can 
be absorbed by x and Co, of which x is a function of time 
(or orbital phase) and Co is unchanged. 

While for general NS-NS or NS-WD binary pulsars, 
77ii i= n%2 and Si ^ 0, Si 7^ 0, the two spins precesses 
at different velocities (Oi 7^ f2 2 ), therefore, S varies both 
in direction and magnitude (S = Si + S 2 ). Then Qq is 
a function of time by Eq(£l} ((Iq ^ 0), which leads to 
the nutation (P(, ^ 0) by EqJJjJ. Therefore, for a general 
NS-NS or NS-WD binary pulsar, there is nutation effect 
beside the precession of the orbit and the apsidal motion. 

Thus, for a general binary pulsar, the three constant 
parameters x obs , ui obs , and P h obs together, can largely 
eliminate the trend of residuals, or the additional time 
delay which can be represented as DM vs time, or DM 
vs orbital phase. 

Since Qq varies with time, then x, Co and P^ also vary 
with time, as shown by EqJTJ, EqJSJ) and EqJSJ respec- 
tively, whereas, the measured x obs , io obs and P b ofcs , are all 
constants, or the average values of the true effects, thus 
for such binary pulsar which have small orbital periods, 

1. e., a few hours (x obs and P b obs vary rapidly), higher or- 
der derivatives, such as x obs and P§ bs are necessary to 
eliminate the trend of residuals. The fitting of the DM 
vs orbital phase in 47 Tuc J in this paper actually in- 
cluded such higher order of derivatives, because x, Co, and 
Pb used in fitting are given by EqlJTJ, EqJSJ) and Eq© 
respectively, which are all functions of time. 

The relationship of the geodetic precession induced sec- 
ular variabilities and the additional Roemer time delay 
as well as the DM variation is summarized in Fig 3. 

Therefore, comparing with the secular variations, x 
and Pb, the variation of DM of 47 Tuc J is just the short- 
term effect (t « Pb/20) of geodetic precession in a binary 
pulsar. 

DM, and secular variabilities, such as x and Pb, have 
been interpreted separately by different models. While 
the geodetic precession provides an unified model which 



6 



can well explain these variabilities (both short term and 
long-term) . 

On the other hand, the DM variations provide new ev- 
idences of the geodetic precession effect in binary pulsar 
systems, beside the secular variabilities, and the theoret- 
ical prediction of the moment of inertia of neutron stars. 

The numerical results of spin angular moment of white 
dwarf star and the pulsar, as well as moment of inertia of 
the pulsar provides very useful information on both the 



structure of white dwarf stars and neutron stars. 

The interpretation of the DM variation by the dynamic 
effects in the three binary pulsars (the frequency prob- 
lems in the previous explanation are solved automati- 
cally) indicts that the structure of the ISM between the 
these pulsars and the Earth might be very stable, which 
is similar as PSR B1855+09 6J. This might provide new 
information in the understanding of ISM. 



TABLE I: Parameters obtained by fitting DM vs time, x obs 
and P£ bs of 47 Tuc J 



XjSi -V/S 2 ^LSi 4>i W i Vi 

0.9192 1.476 0.9203 0.1095 3.594 0.3640 

£i(g cm z s _i ) S2(g cm z s _i ) mi rn~2 e i 

5.985 x 10 48 4.933 x 10 49 1.949 0.0264 1.929 x 10" 5 1.570 



By Si the moment of inertia of the pulsar can be obtained, 

ii = 2.0 x 10 45 g cm 2 . Notice that mi and rri2 are in sun 
mass, all angles are in radian. The errors of the values in the 
above table is about 30% 
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FIG. 1: Binary geometry and definitions of angles. The in- 
variable plane (x-y) is perpendicular to the total angular mo- 
mentum, J. The inclination of the orbital plane with respect 
to the invariable plane is which is also the precession 

cone angle of L around J. The orbital inclination with re- 
spect to the line of sight is i. r/o is the phase of the orbital 
plane precession, lo is longitude of periastron (point P). 
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FIG. 2: The measured@ and predicted DM as a function of 
orbital phase in 47 Tuc J . The measured upper plot represents 
the measurement made with 1999-2002 high-resolution data, 
displayed O.Olcm" 3 pc below their measured values for clarity. 
The second plot represent the DM obtained using only the 
data from the best observation (1999 October 11), displayed 
0.02cm -3 pc below their measured values0- The solid curve 
of the upper plot represent the best fit. And the solid curve 
of the second plot is obtained by setting the time 0.3714 year 
later than the upper fitting. 




FIG. 3: Geodetic precession leads to an additional orbital 
precession velocity, Oo, which can be absorbed by three pa- 
rameters, x, Co and P b . Which in turn cause an additional 
time delay, A^, which can be attributed to DM variation. 



